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>D ■ Lattice-based investigations of two fundamental QCD quantities are described, 

namely the gluon and quark propagators in Landau gauge. We have studied the 

Cn ' Landau gauge gluon propagator using a variety of lattices with spacings from 

^ a = 0.17 to 0.41 fm to explore finite volume and discretization effects. We also 

Q\ , introduce the general method of "tree-level correction" to minimize the effect of lat- 

£*\J ■ tice artefacts at large momenta. We demonstrate that it is possible to obtain scaling 

f"^ ■ behavior over a very wide range of momenta and lattice spacings and to explore 

^^ ' the infinite volume and continuum limits of the Landau-gauge gluon propagator. 

f**^ ' These results confirm the earlier conclusion that the Landau gauge gluon propa- 
gator is infrared finite. We study the Landau gauge quark propagator in quenched 

/~~s QCD using two forms of the C(a)-improved propagator with the Sheikholeslami- 

-«^ i Wohlert quark action with the nonperturbative value for the clover coefficient c sw 

^ C"| and mean-field improvement coefficients in our improved quark propagators. We 

O i' again implement an appropriate form of tree-level correction. We find good agree- 

I ment between our improved quark propagators. The extracted value of infrared 

C-1h quark mass in the chiral limit is found to be 300 ± 30 MeV. We conclude that 

Qj i the momentum regime where the transition from nonperturbative to perturbative 

^- QCD occurs is Q 2 ~4 GeV 2 . 
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1 Introduction 

Lattice gauge theory is currently the only known "first principles" approach to 
studying nonperturbative QCD. It is therefore important for lattice QCD ta 
provide constraints and guidance for the construction of quark-based modelsJ 
and to provide an indication of the momentum regime at which we can expect 
perturbative QCD to become applicable. The quark and gluon propagators are 
two of the most fundamental quantities in QCD. There has been considerable 
interest in the infrared behavior of the gluon propagator as a probe into the 
mechanism of confinement and by studying the scalar part of the quark prop- 
agator, the mass function, we can gain insight into the mechanisms of chiral 
symmetry breaking. Both are used as input for other quark- model calculations. 





Dimensions 


P 


a (fm) 


Volume (fm 4 ) 


Configurations 


lw 


16 3 x 32 


5.70 


0.179 


2.87 3 x 5.73 


100 


li 


16 3 x 32 


4.38 


0.166 


2.64 3 x 5.28 


100 


2 


10 3 x 20 


3.92 


0.353 


3.53 3 x 7.06 


100 


3 


8 3 x 16 


3.75 


0.413 


3.30 3 x 6.60 


100 


4 


16 3 x 32 


3.92 


0.353 


5.65 3 x 11.30 


100 


5 


12 3 x 24 


4.10 


0.270 


3.24 3 x 6.48 


100 


G 


32 3 x 64 


6.00 


0.099 


3.18 3 x 6.34 


75 



Table 1: Details of the lattices used to calculate the gluon propagator. Lattices lw and li 
have the same dimensions and approximately the same lattice spacing, but were generated 
with the Wilson and improved actions respectively. Lattice 6 was generated with the Wilson 
action. 



2 Gluon Propagator 

We use an 0(a 2 ) tree-level, tadpole-improved actiorH and for the. tadpole 
(mean-field) improvement parameter we use the plaquette measurda. A full 
description and discussJmi_of the gluon propagator results summarized here 
can be found elsewhereoEm 

Gauge fixing on the lattice is achieved by maximizing a functional, the 
extremum of which implies the gauge fixing condition. The usual Landau 
gauge fixing functional implies that ^ d^A^ = up to 0(a 2 ). To ensure 
that gauge dependent quantities are also O{o?) improved, we implement the 
analogous 0{a 2 ) improved gauge fixingQ The dimensionless lattice gluon field 
Afj,(x) is calculated from the link variables in the usual way, which agrees with 
the continuum to 0(a 2 ). We then calculate the scalar part of the propagator 



D(x-y) = J2<Mv)M*))- (!) 



To isolate the nonperturbative behavior of the gluon propagator, we can divide 
the propagator by its lattice tree level form (i.e., that of lattice perturbation 
theory) U For the momentum space gluon propagator D(q 2 ), we see that in the 
continuum q 2 D(q 2 ) will approach a constant up to logarithmic corrections as 
q 2 — > oo because of asymptotic freedom. The continuum tree-level propagator 
is l/q 2 . We also expect asymptotic freedom on the lattice despite finite lattice 
spacing artefacts. We define the lattice q^ such that the lattice D tT00 (q) = l/q 2 , 
and use this momentum throughout. This is referred to as tree- level correction 
and we have seen that it significantly reduces discretization arrors at large 
momenta. For the two actions considered here, this means that we work with 



the momentum variables defined as 



-sin— , <^-V Sln (— j + 3 Sm (— J' (2) 

for the Wilson and improved actions respectively. All figures (quark and gluon 
propagators) have a cylinder cut imposed upon them, i.e. all momenta must 
lie close to the lattice diagonal. In Table [l| we show the various lattices that 
we have studied for the gluon propagator. In Fig. [y we plot q 2 D(q 2 ) for a fine 
unimproved Wilson action and for our finest improved action. Despite having 
very different lattice spacings the agreement is excellent for the entire interme- 
diate and high-momentum regime. The small discrepancy in the deep infrared 
due to finite volume effects is not apparent in this way of plotting that data. 
We plot D(q 2 ) for five different lattice in Fig. g and see pleasing agreement for 
the results. Note that we are plotting bare quantities only and there is thus an 
overall wavefunction renormalization for the gluon propagator (i.e., Z3(/i,a) 
for the renomalization point fi). The vertical scale is thus unimportant and 
only the variation with momentum is relevant. This way of presenting the 
data shows that there is a small residual finite volume dependence, where the 
infrared gluon propagator is decreasing with increasing lattice volume V. We 
have performed a fit as a function of 1/V and have seen that the large vol- 
ume P = 3.92 lattice gives results which are very close to the infinite volume 
limit. In Fig. |we plot D(q 2 ) in the intermediate and ultraviolet regime and 
compare with the three-loop perturbative QCD form. We see that the above 
q ~ 2 GeV the agreement is excellent, but that below this momentum scale 
nonperturbative effects are becoming apparent. 



3 Quark Propagator 

The Landau gauge quark propagator results summarized here have been pre- 
sented and discussed in more detail elsewhereEj All 0(a) errors in thcicrmion 
action can be removed by adding appropriate terms to the Lagrangianxa. It is 
then usual to perform appropriate field transformations to improve the quark 
operators as wellla 

In the continuum, the quark propagator has the following general form, 

stn) = - = zc{p) (n 

[P) i M c (p) + B C { P ) i i> + M c (p) ' l ' 

We expect the lattice quark propagator to have a similar form, but with ]k 
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Figure 1: Comparison of the gluon propagator from the finest improved lattice (lattice 
li, /9 = 4.38) and the finest Wilson lattice (lattice 6, /3 = 6.0). Data has been cylin- 
der cut and the appropriate tree-level corrections have been applied. The data from lat- 
tice 6 is half-cut whereas lattice li displays the full Brillouin zone. We have determined 
Z3(improved)/Z3(Wilson) = 1.08 by matching the vertical scales of the data. 



replacing p": 



S(p) = 



z(p) 



i ft + M(p) 



(4) 



where k is a new 'lattice momentum', k^ = i sin(p M a). We do not have 
sufficient space here to describe the hybrid tree-level correction that was used 
for the quark propagator., results presented here, but a detailed description 
has recently been givenlij We again use the cylinder cut to further reduce 
hypercubic discretization artefacts. As for the gluon propagator the results 
for Z(p) are for the bare quantity only and contain an overall renormalization 
constant Zi{\i,ci). In Fig. [| the vertical scales have been adjusted so that 
the two sets of results are renormalized and hence coincide at 2.1 GeV. In 
this figure we see the charactreistic dip in the infrared for Z(p), which occurs 
also in model Dyson-Schwinger equation studiefJ of dynamical chiral symmetry 
breaking. This dip has essentially disappeared by around 2 GeV. The improved 
action correspondning to Sr is the one we prefer and it gives the more expected 
ultraviolet behavior of the renormalized Z(p), i.e., it tends towards a constant. 
In Fig. @ we see the characteristic behavior of the quark mass function familiar 
from quark model studies! and that the transition to the perturbative regime 
is occuring at approximately 2 GeV. Note that there is no renormalization of 
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Figure 2: Comparison of the gluon propagator generated with an improved action on five 
different lattices. We find good agreement down to q ~ 500 MeV. At the lowest accessible 
momenta the data points drop monotonically with increasing volume, but the lowest point 
(on the largest lattice) shows signs of having converged to its infinite volume value. For 
comparison with perturbation theory, a plot of the continuum, tree-level gluon propagator 
(i.e., 1/q 2 appropriately scaled) has been included. 



the quark mass function. At large momenta the mass function should become 
the running quark mass of perturbative QCD. Finally in Fig. ra we present a 
simple quadratic extrapolation to the chiral limit for the available Si data. The 
slight dip is not statistically significant and is almost certainly a residual lattice 
artefact. The infrared mass (i.e., at p = 0) in the chiral limit is approximately 
330 ± 30 MeV, which is characteristic of the constituent quark mass scale. 



4 Summary and Conclusions 

The gluon propagator has been calculated on fine unimproved lattices and 
on a variety of improved lattices with an 0(a 2 ) improved action in 0(a 2 ) 
improved Landau gauge. The infrared behavior of this propagator strongly 
suggests the the Landau gauge gluon propagator is infrared finite. We have 
ruled out the l/<? 4 behavior popular in some Dyson-Schwinger quark model 
studies and indeed any infrared singularity appears to be very unlikely. The 
possible effects of lattice Gribov copies remains a very interesting question 
and we are currently carrying out similar studies across a variety of lattices 
in Laplacian gauge, which is a Landau-like smooth gauge fixing, but is free of 
Gribov copies. 
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Figure 3: Comparison of the lattice gluon propagator with that obtained from perturbation 
theory, in the ultraviolet to intermediate regime. The continuum expressions are tree-level 
(i.e., 1/q 2 appropriately scaled) and the three-loop perturbative QCD expression. 



We have used two different definitions of the 0(a) improved quark propa- 
gator, corresponding to the quark propagators denoted Si and Sr. We make 
use of asymptotic freedom to factor out the tree level behaviour, replacing it 
with the 'continuum' tree level behaviour Z(p) = 1, M(p) = m. This tree-level 
correction dramatically improves the data. We find that M(0) approaches a 
value of 300 ± 30 MeV in the chiral limit, which is very much in keeping with 
the concept of a "constituent quark mass" and agrees with the infrared values 
of the quark mass commonly used in model studieslil We also find a significant 
dip in the value for Z(p) at low momenta. This is again entirely consistent 
.with what is found in model studies of dynamical chiral symmetry breaking 
El. An examination of Figs. 01 0, 0, H, and provide a clear indication that 
perturbative QCD behavior is not becoming dominant in the gluon and quark 
propagators until we reach momenta of order Q 2 ~ 4 GeV 2 . 
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Figure 4: The tree-level corrected Z(p) for c sw = NP and k = 0.137. The figure shows Z(p) 
vs. momentum in physical units, and after rescaling ( "renormalizing" ) so that Z(2.1GeV) = 
1. The infrared agreement after rescaling is very good and we take the tree-level corrected 
Z(p) from Sr with the nonperturbative c sw as the best estimate for this quantity. 



9. 
10. 



D.B. Leinweber, J.I. Skullerud, A.G. Williams, and C. Parrinello, Phys. 

Rev. D 58, 031501 (1998); ibid. 60, 094507 (1999); Erratum - ibid. 61, 

079901 (2000); D.B. Leinweber, C. Parrinello, J.I. Skullerud, and A.G. 

Williams Nucl. Phys. B (Proc. Suppl.) 73, 629-631 (1999); ibid., 626-628 

(1999). 

F.D.R. Bonnet, P.O. Bowman, D.B. Leinweber, and A.G. Williams, 

Phys. Rev. D 62, 051501 (2000); ibid. 64, 034501 (2001). 

M. Liischer et al, Nucl. Phys. B 478, 365 (1996). 

C. Dawson et al, Nucl. Phys. B (Proc. Suppl.) 63, 877 (1998). 

G. Heatlie et al, Nucl. Phys. B 352, 266 (1991). 

J.I. Skullerud and A.G. Williams, Nucl. Phys. B (Proc. Suppl.) 83- 

84, 209-211 (2000); Phys. Rev. D 63, 054508 (2001); J.I. Skullerud, 

D.B. Leinweber, and A.G. Williams, to appear in Phys. Rev. D, h.ep- 



lat/0102013 



-^- 



c m = 91 MeV 
x m = 54 MeV" 
o m = 35 MeV 



*o ; * Xx » ; > ; * :<xx *> <xxxxx»<xxxxx>> 



2 4 

P (GeV) 





,_ 




1 


1 ' 


1 , i i i | i i , i 




ft 








Sr. c sw =np : 


400 




I 






- 


300 


- 








- 


200 
100 


7 




4 


s 

4 


o m=91 MeV - 

** 



2 4 

P (GeV) 



Figure 5: The tree-level corrected M(p), for c sw = NP, using Sj (left) and Sr (right). 
We find good agreement between the two data sets, both in the infrared and ultraviolet. 
The residual disagreement at intermediate momenta is a pointer to lattice artifacts that we 
have not brought under full control, even with nonpcrturbative improvement and tree-level 
correction. 
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Figure 6: The tree-level corrected mass function from Sj with c sw = NP, with the bare mass 
extrapolated to zero using a quadratic fit. The small dip at p ~ 1.6 GeV is not statistically 
significant and may be due to residual lattice artifacts. The non-zero values for M(p) in 
the chiral limit are entirely due to dynamical chiral symmetry breaking and provide a direct 
measure of this effect. 



